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We study the medium modification of jets correlated with large transverse momentum photons at 
the LHC via a Monte-Carlo transport model which incorporates both elastic collisions and radiative 
energy loss experienced by the parton showers. Numerical results are presented for the modification 
of tagged jet yield, the photon-jet momentum imbalance, and their azimuth separation distribution. 
The use of photon-jet correlations as jet quenching tomography tools is addressed by studying 
the photon-jet events with different values of xt = Pt,j/pt,t It is found that larger xt jets are 
mainly produced from the surface of the medium, while smaller xt jets are from those jets which 
have traversed longer distance of medium. Combined with the momentum spectrum of the photon- 
tagged jets, this leads to stronger suppression and centrality dependence for larger xt jets, and 
weaker suppression (or enhancement) for smaller xt jets. 



The energetic probes produced in early hard scattering 
processes provide valuable tools for studying the prop- 
erties of the quark-gluon plasma (QGP) created in rcla- 
tivistic heavy-ion collisions. Large transverse momentum 
(pt) jets arc particularly useful as they directly inter- 
act with the constituents of the dense QCD matter that 
they traversed, and are modified/quenched in the pro- 
cess Jet quenching has been confirmed by multiple 
experimental observations, such as the significant sup- 
pression of single inclusive high px hadron production 
in nucleus-nucleus collisions compared to binary collision 
scaled proton-proton collisions at the same energies 0-11] . 

A lot of effort has been devoted to understand the mod- 
ification of jets by the hot and dense strongly-interacting 
QCD matter. Jet modification originates from a com- 
bination of drag, diffusion and stimulated gluon emis- 
sion experienced by the partonic jets while propagat- 
ing through the highly excited nuclear medium [5|-l8|. 
Various transport coefficients, e.g., e = dE/dt, 62 = 
d{AE)'^/dt and q = d{ApT)'^ /dt, have been employed 
to quantify these modification effects d, There are 
currently a number of perturbative QCD based schemes 
for the treatment of the radiative part of parton energy 
loss [lol - [l^ . A systematic comparison of these differ- 
ent jet quenching schemes has been carried out in Ref. 

in terms of the "brick" problem. Many phenomeno- 
logical calculations have been performed to describe the 
observed jet modification measurements, such as the sup- 
pression of single inclusive high px hadron production 
[16l - [20| , and the nuclear modification of dihadron [2l|, [2^ 
and photon-hadron p3l - [25| correlations in high energy 
nucleus-nucleus collisions. 

Recently, sophisticated experimental techniques have 
been developed to reconstruct the full jets emitted in rela- 
tivistic heavy-ion collisions. This provides more stringent 
test on our understanding of jet modification in terms of 
both leading and subleading fragments of the jet showers 



[26|. Various Monte-Carlo generators are being devel- 
oped to simulate the modification of jet propagation in 
dense nuclear medium [27, 30] . The study of in- medium 
evolution of the full jets helps to understand how energy 
and momentum arc deposited by the jet showers and how 
the medium responds to jet propagation and energy de- 
position j3ll - l33| . 

Meanwhile, the much larger kinematics achieved at 
the LHC now allows us to investigate medium effects on 
jets with very large transverse energies (over a hundred 
GeV). One of the first exciting full jet results from the 
LHC is the strong modification of the momentum im- 
balance distribution between the correlated jet pairs in 
central Pb-|-Pb collisions, while no strong modification 
is observed for the dijct relative azimuth separation dis- 
tribution [H, [3^. This result indicates significant en- 
ergy loss experienced by the subleading jets during their 
propagating through the hot and dense medium created 
in Pb+Pb collisions. There have been various models 
employed to explain the observed dijets asymmetry in 
Pb-hPb collisions at the LHC 



Another promising observable is the modification of 
the jets correlated with photons or electrowcak bosons 
in relativistic heavy-ion collisions. Such channel has 
been regarded as the "golden" channel for jet quench- 
ing study due to the fact that the neutral boson trig- 
gers, once produced, will escape the medium without 
further interaction. Thus a stronger constraint is pro- 
vided by the triggers for the away-side tagged jets that 
we are interested in [i^. Furthermore, many biases 
that are present in single inclusive jet/hadron suppres- 
sion and dijet/dihadron/jet-hadron correlations are im- 
proved, such as the deep falling spectra, trigger bias, and 
the sensitivity to the background fiuctuations. Experi- 
mental measurements of photon-triggered fragmentation 
function at RHIC [H, and the momentum imbalance 
of correlated photon-jet pairs at the LHC [i^ are consis- 
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tent with parton energy loss calculations [23l - l25l . 1471 . |48| . 

In this work, we focus on studying the medium mod- 
ification of jets correlated with high pT photons at the 
LHC energies. A Monte-Carlo transport code is devel- 
oped to simulate in-medium jet shower evolution and jet 
modification in dense QCD medium, with the inclusion 
of both stimulated gluon radiation and the elastic colli- 
sions experienced by the shower partons. The effect of 
elastic collisions is quantified by a few transport coeffi- 
cients, namely the longitudinal drag and diffusion as well 
as transverse momentum broadening, while the radiative 
energy loss is simulated by employing the single gluon 
emission spectrum obtained from Higher Twist jet en- 
ergy loss calculation [3 ■ In section II, we provide 
some details of how the simulation of jet shower trans- 
port in medium is performed. The numerical results for 
photon-jet correlations in Pb+Pb collisions at the LHC 
energies are presented in section HI, where we address the 
use of photon-jet correlations as jet quenching tomogra- 
phy tools by studying the tagged jets with different xt 
values. The last section contains our summary. 



I. IN-MEDIUM JET TRANSPORT 

In this section, we describe our Monte-Carlo method 
for simulating the transport and modification of a par- 
tonic jet shower in dense QCD matter. We assume each 
parton of the jet shower follows classical trajectories be- 
tween each jet-medium interaction that induces the mo- 
mentum change of the parton (either by an elastic colli- 
sion or the radiation of a gluon) . We keep track of all the 
propagating partons and the radiative gluons until they 
traverse the dense medium. 

The probability for a parton at each time step Ai to 
have an momentum exchange by an elastic collision with 
the medium constituents is determined by. 



Pcon(i,At) = At/t 



mf 



(1) 



where imf is an "effective" mean free time, thus f/tmf 
is the rate of the elastic collisions (energy and momen- 
tum exchange). If there is a jet-medium interaction 
through an elastic collision, the energy and momentum 
exchange is sampled from a probability distribution of 
P(Ai;,Apj.,i,At). 

Generally, the above mean free time trai and the dis- 
tribution of the momentum and energy exchange can 
be determined if the details about the medium struc- 
ture and jet-medium interaction are known. In this ap- 
plication, we assume that the probability distribution 
P{AE, Ap±,t, At) is a Gaussian distribution in both lon- 
gitudinal and transverse directions, with the mean and 
invariance provided by the longitudinal drag etmf and 
diffusion e2imf , and transverse broadening qtmi- 

The Gaussian form for the momentum exchange 
through the elastic collisions originates from the assump- 
tion of multiple scatterings with small momentum trans- 
fer for each scattering. Such limit corresponds to setting 



the mean free time equal to the time step (imf = At), i.e., 
the Langevin approach. The variation of t^[ values does 
not affect the results for Gaussian form of momentum ex- 
change and large extent of the medium {At < imf ^ L). 
Realistic values of mean free time <mf may be used in the 
future study when the method is extended to include the 
case when the momentum exchange is a non-Gaussian 
form, such as few hard collisions with a power-law tail. 

We simulate the radiation of gluons by employing the 
single gluon emission spectrum from the higher twist cal- 
culation [Ulil, 



dxdli dt 
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where x and l± arc the fractional energy and transverse 
momentum carried by the radiated gluon. P{x) is the 
vacuum splitting function, and tform = 2Ex{l — x)//^ is 
the gluon formation time. 

The probability for a parton in each time step At to 
radiate a gluon is determined by the average number of 
gluons: 



Prad(i, At) = {Ng{t,At)) = At / dxdl 



dN„ 



dxdl\dt 



(3) 



If there is a gluon emitted in a given time step, the energy 
and the momentum of the radiated gluon are sampled 
from the following probability distribution. 



P(x,?i,t, At) 



At dNg 
(Ng) dxdl\dt 



(4) 



To simulate the medium modification of jet shower, 
we need to provide the initial conditions before the jet- 
medium interaction. In this application, the production 
of the photon-jet pairs from early hard scattering pro- 
cesses are generated from PYTHIAS. I [i^ for p-|-p col- 
lisions at ^JSNN = 2.76TeV. The full jet reconstruction 
is performed via FASTJET package [50| from which one 
can find the detailed content of the full jets, such as the 
energy and momentum for each parton in a jet. Such 
initial conditions are provided to our Monte-Carlo code 
for the further simulation of in-medium transport and 
medium modification. 

The initial photon-jet production points is sampled ac- 
cording to the distribution of the binary collisions sim- 
ulated from the Glauber Model. The evolution pro- 
files (energy /entropy density, temperature and fiow ve- 
locities) of the bulk QGP medium that jets interact 
with is provided by hydrodynamic simulation. In this 
work, we employ the (2+I)D viscous hydrodynamics 
model (VISH2+ij developed by The Ohio State Univer- 
sity group j5ll - l53{ . with two-component Glauber model 
for hydrodynamics initial conditions. The code version 
and parameter tunings for Pb-|-Pb collisions at LHC en- 
ergies are taken as in Ref. [53| . 

When simulating the propagation of the jet shower 
through the hydrodynamic medium, we only allow the 
vacuum radiated gluons to interact with the medium 
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when their formation times tform are reached. We impose 
a minimum energy cutoff ttT for the medium-induced ra- 
diation spectrum in order to take into account the bal- 
ance between gluon radiation and absorption from the 
hydrodynamic medium. When the local temperature 
of the medium drops below the transition temperature 
160MeV, jets are assumed to decouple from the medium. 
After jets traverse the hydrodynamic medium, we provide 
the detailed information of the modified jets to FAST- 
JET package to perform the jet reconstruction and then 
compare with the vacuum jets. In this work, we have not 
included the hadronization process which we postpone to 
an upcoming effort. 

As mentioned above, the jet shower evolution in dense 
QGP medium are controlled by three transport coeffi- 
cients e, 62 and q. Thse coefficients can be calculated 
if the detailed information about medium structure and 
jet- medium interaction is known [s^ ]. Here we assume 
that the diffusions in transverse and longitudinal direc- 
tions are symmetric, and the fluctuation-dissipation re- 
lation may be applied between the drag and diffusion, 
q » 2e2 ~ 4Te. The functional form of the transport 
coefficients is taken as g oc T"^, with the overall normal- 
ization to be fitted to one data point. The sensitivity 
of jet in-medium transport and modification to different 
parameterizations of jet transport coefficients will be in- 
vestigated in the future study. 



II. NUMERICAL RESULTS FOR PHOTON-JET 
CORRELATIONS 

In this section, we present some numerical results for 
the medium modification of jets correlated with high 
Pt photons in Pb-|-Pb collisions at ^/sWn = 2.76TeV 
at the LHC. One of the most interesting quantities is 
the distribution of the tagged jet momentum fraction 
XT ~ Pt.j /PT,-y, which may be defined as 



f{xT) = {l/N^)dNj^/dxT 

P{xt) = {l/Nj^)dNj^/dxT 



(5) 



The above two definitions differ only by the normaliza- 
tion. Given the kinematic cuts for the trigger photons 
and associated jets, the former distribution is normalized 
to the fraction of photons having the associated jet pair 
Rjj and the latter gives unity when integrating out the 
momentum fraction xt- The nuclear modification factor 
Iaa for the tagged jet distribution is usually defined as 



Iaa{xt) = fAA{xT)/fpp{xT) 



(6) 



In Fig. [T] we show the zx-integrated nuclear modi- 
fication factor Iaa — R^^/R^f^ (left) and the average 
fractional momentum loss (Axt) = {xt)\pp — {xt)\aa 
(right) as a function of centrality (the participant num- 
ber A'part)- The data points shown here are obtained 
from CMS measurements of Rj^ and (xt) for Pb-|-Pb 
coUisions and PYTHIA-f HYDJET references In 
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FIG. 1: (Color online) The nuclear modification factor Iaa 
and the average momentum loss {Axt) for photon-triggered 
jets as a function of centrality for Pb-fPb collisions at the 
LHC. The jet size is i? = 0.3. 



this work, the CMS kinematic cuts are taken for all the 
presented results: the photon momentum pT,7 > 60GeV, 
jet momentum px.j > 30GeV, and their azimuth sepa- 
ration A(j> = |(/)-, — (pjl < Itt/S. The size of the jets is 
taken as i? = 0.3 in this plot. The model parameter q 
(accordingly e and 62) is fitted to the value of Iaa ~ 0.7 
in most central (0 — 10%) Pb-|-Pb collisions; we obtain 
the gluon transport coefficient go = 6.5GeV^/fm at hy- 
drodynamics initial time tq = 0.6fm/c. Due to stronger 
medium effect from peripheral to central collisions, one 
observes the decrease of the nuclear modification factor 
Iaa and the increase of the average fractional momentum 
loss (Axt) for the photon-triggered jets. 
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FIG. 2: (Color online) The distribution of the momentum im- 
balance variable xt between triggered photons and associated 
jets for most central (0 — 10%) Pb-f Pb collisions at the LHC. 
The jet size is R = 0.3. 

One may look at the distributions of the momentum 
imbalance variable xj between triggered photons and the 
away-side jets. In Fig. [51 the normalized distribution 
P{xt) is shown for both p+p and most central (0 — 10%) 
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FIG. 3: (Color online) The distribution of the relative az- 
imuth separation Atj) distribution between triggered photons 
and associated jets for most central (0 — 10%) Pb+Pb colli- 
sions at the LHC. The jet size is R = 0.3. 



Pb-I-Pb collisions at the LHC. One observes the shift of 
the distribution to smaller xj values for the distribution 
in Pb-fPb collisions compared to p+p collisions. This 
indicates that the tagged jets lose energy during its prop- 
agation through the hot and dense QGP created in the 
Pb-|-Pb collisions. 

The azimuth separation between the trigger photons 
and associated jets Acj) = \(j)-y ~ <j>j \ is another interesting 
quantity. In Fig. [3l the normalized distribution P{A(f)) is 
shown for both p+p collisions and most central (0 — 10%) 
Pb-I-Pb collisions. While there is strong modification for 
the momentum imbalance xj distribution as shown in 
Fig. [21 no strong modification is observed for photon-jet 
azimuth separation distribution. This is partly due 
to the kinematic cuts applied on the selection of photon- 
jet events. Jets which are more strongly deflected by the 
dense medium are those having experienced larger en- 
ergy loss. Thus after applying the kinematic cuts, many 
of such associated jets are not selected in the final re- 
sult. Such effect is illustrated by the dashed curve which 
represents no px cut applied to the associated jet mo- 
mentum after the in-medium evolution. One observes 
some broadening on the photon-jet azimuth separation 
A(j) distribution. Such broadening effect is further cut off 
by the constraint on the photon-jet azimuth separation 
Acj) > 7 /Sir. 

In Fig. m the photon-triggered nuclear modification 
factor Iaa is plotted as a function of the momentum 
fraction xt for Pb+Pb collisions at the LHC. Results 
for different centralities are shown for comparison. One 
may observe that Iaa decreases from lower xt values to 
higher xt values. This means that the yield for the asso- 
ciated jets with high xt values is suppressed in Pb+Pb 
collisions due to the interaction with the QGP medium, 
while the yield for lower xt jets is enhanced due to the 
fact that the momentum fraction variable xt distribution 
shifts from higher xt to lower xt values. We observe that 
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FIG. 4: (Color online) The nuclear modification factor Iaa for 
photon-triggered jets as a function of xt for Pb+Pb collisions 
at the LHC in different centralities. The jet size is R — 0.3. 
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FIG. 5: (Color online) The nuclear modification factor Iaa 
for photon-triggered jets as a function of xt for most central 
(0 - 10%) Pb+Pb coUisions at the LHC. Results for different 
jet sizes are compared. 



such medium modification effect increases from periph- 
eral collisions to central collisions. 

The above modification effect on the photon-triggered 
jets depends on the jet size. In Fig. [5] the photon- 
triggered nuclear modification factor Iaa is plotted as 
a function of xt for most central (0 — 10%) Pb+Pb col- 
lisions at the LHC. Results for different jet sizes {R = 
0.2,0.3,0.4,0.5) are shown for comparison. We observe 
that the suppression for the associated jet yield is similar 
for high XT values for different jet sizes, and for smaller 
XT values the enhancement is larger for larger jet size. 

In the above, we have presented the medium modifi- 
cation of photon-triggered jet yield, the photon-jet mo- 
mentum imbalance, and their azimuth separation distri- 
bution for Pb+Pb collisions at the LHC. Now we ad- 
dress the use of photon-triggered jets as a tomographic 
tool for studying jet quenching in relativistic heavy-ion 
collisions. In particular, we may classify the photon-jet 
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FIG. 6: (Color online) The density distribution of the initial 
photon-jet production points (xinijj/ini) when the triggered 
photons are taken along out-of-plane directions {\(f>-^ — n /2\ < 
it/12) for Pb+Pb collisions at the LHC. Jets with different 
final XT values are compared: xt = [0.5,0.6] for (a) and (c); 
XT = [0.9, 1] for (b) and (d). Two different centralities are 
compared: - 10% for (a) and (b); 20 - 30% for (c) and (d). 
The jet size is _R = 0.3. 



events according to different values of final xt . This al- 
lows us to use photon-triggered jets to probe different 
regions of the hot and dense QGP created in the col- 
lisions. To illustrate this, we show in Fig. [6] the den- 
sity distribution of initial photon-jet production points 
(a^iniij/ini) in thc trausvcrsc plane. Here wc take those 
events with thc triggered photons that propagate along 
out-of-planc directions — 7r/2| < 7r/12). Wc compare 
the initial photon-jet production point distribution for 
jets with different xt values: xt ~ [0.5, 0.6] for (a) and 
(c); XT = [0.9,1] for (b) and (d). Results for two differ- 
ent centralities are also shown for comparison: — 10% 
for (a) and (b); 20 - 30% for (c) and (d). The jet size is 
taken ds, R = 0.3 in this plot. 

If one takes photon-jet events for all values of xt and 
propagating directions, the density distribution of initial 
photon-jet production points would follow binary colli- 
sion distribution (if no kinematic cuts arc applied after 
medium evolution). When selecting jet events with cer- 
tain values of xt , one can see that those jets with larger 
values of final xt are mainly produced at the surface of 
the medium thus have traversed short medium length, 
while smaller xt jets are more from these surviving jets 
which have traversed much longer distance. Combined 
with selecting trigger photons with different directions, 
we may probe different regions of the hot and dense 



medium. Thus photon-triggered jets with different xt 
values provide very useful tomographic tools for study- 
ing jet modification in relativistic heavy-ion collisions. 
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FIG. 7: (Color online) The average values of the jet initial 
production points (j/ini) when the triggered photons are taken 
along out-of-plane directions {\(t)^ — n/2\ < n/12) for Pb-|-Pb 
collisions at the LHC. Results for different centralities are 
compared. The jet size is R = 0.3. 
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FIG. 8: (Color online) The average values of the jet initial 
production points (j/ini) when the triggered photons are taken 
along out-of-plane directions {\cj)~/ — tt/2\ < 7r/12) for central 
- 10% Pb-I-Pb coUisions at the LHC. Results for different 
jet sizes are compared. 

The above effect can be made more quantitative. In 
Fig. [3 we show the average values for the photon-jet 
initial production points (j/ini) ff one looks at the trig- 
gered photons that propagate along out-of-planc direc- 
tions (10^ — 7r/2| < 7r/12). One may see that on average 
jets with smaller xt values may have traversed about a 
few fm longer distance than larger xt jets. It is noted 
that such extra traveled path is the hottest part of the 
QGP medium. We observe that such effect is larger for 
central collisions than peripheral collisions. 

The jet size dependence of the above effect is shown 
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in Fig. [U where the average values of the jet initial pro- 
duction points (j/ini) is plotted as a function of xt for 
most central (0 - 10%) Pb+Pb collisions at the LHC. 
The triggered photons is taken to be along out-of-plane 
directions — tt/2\ < 7r/12) as above. One observes 
that the above effect show weak dependence on the jet 
size (stronger for larger jet size). 
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FIG. 9: (Color online) The nuclear modification factor Iaa 
for photon-triggered jets as a function of centrality for Pb-f Pb 
collisions at the LHC. Results for different xt values are com- 
pared. The jet size is R = 0.3. 

Since the tagged jets with different final Xt values are 
mainly from different regions of QGP and traverse dif- 
ferent lengths of the medium, they show different sensi- 
tivity to parton energy loss. In Fig. 1^1 we show nuclear 
modification factor Iaa for the photon-triggered jets as a 
function of centrality (iVpart) for Pb-|-Pb collisions at the 
LHC. Results for different values of xx bins are shown 
for comparison. The jet size is taken as i? = 0.3 in 
this plot. Due to the combinational effects of different 
traveling distances and the tagged jet spectra, one ob- 
serves different centrality dependence for the modifica- 
tion of photon-triggered jets with different values of xt- 
Jets with larger values of xt show larger suppression and 
stronger centrality dependence, while smaller suppression 
(or enhancement) is observed for smaller xt jets. 



experienced by the jet shower partons is encoded by a 
few non-perturbative transport coefficients such as the 
drag and diffusion, while the radiative energy loss is sim- 
ulated by employing the single gluon emission spectrum 
from Higher Twist jet energy loss formalism. 

With our event-by-event jet shower evolution code, 
we have presented the numerical results for the nuclear 
modification of photon-triggered jet production yield, 
the photon-jet momentum imbalance and their azimuth 
separation distribution. We have addressed the poten- 
tial of photon-jet correlations as a tomographic tool 
for studying jet quenching in heavy-ion collisions by 
classifying photon-tagged jets with different values of 
the momentum fraction xt- We found that jets with 
larger xt values are mainly produced from the surface 
of the medium, while jets with smaller xt values jets 
are mainly from these jets which have traversed longer 
distance. This makes photon-triggered jets very useful 
tomographic tools for studying jet modification in dense 
matter created in relativistic nuclear collisions. Com- 
bined with the tagged jet momentum spectrum, we found 
different sensitivity to parton energy loss for different 
photon-tagged jets: one observes stronger suppression 
and centrality dependence for larger xt jets, and weaker 
suppression (or enhancement) for smaller xt jets. 

The above analysis may be directly applied for study- 
ing the nuclear modification of high px dijets and single 
inclusive jets. The inclusion of the hadronization pro- 
cess the for final partonic jet showers will enable us to 
directly compare to the experimental measurements of 
jet fragment profiles. The exploration of the sensitivity 
of jet modification to the properties of jet transport co- 
efficients will help understand the detailed structure and 
properties of the dense QGP medium. The energy and 
momentum deposition profiles obtained from the simu- 
lation of in-medium jet shower evolution may be com- 
bined with hydrodynamics simulation to investigate the 
response of the medium to jet transport and the influence 
of hydrodynamic background and its fluctuations on full 
jet observables. These will be left in the follow-up effort. 
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